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ABSTRACT 

We present new multicolor photo-polarimetry of stars behind the Southern 
Coalsack. Analyzed together with multiband polarization data from the litera- 
ture, probing the Chamaeleon I, Musca, p Opiuchus, R CrA and Taurus clouds, 
we show that the wavelength of maximum polarization (Xmax) is linearly corre- 
lated with the radiation environment of the grains. Using Far-Infrared emission 
data, we show that the large scatter seen in previous studies of Xmax as a function 
of Ay is primarily due to line of sight effects causing some Ay measurements to 
not be a good tracer of the extinction (radiation field strength) seen by the grains 
being probed. 

The derived slopes in Xmax vs. Ay, for the individual clouds, are consistent 
with a common value, while the zero intercepts scale with the average values of 
the ratios of total-to-selective extinction (Ry) for the individual clouds. Within 
each cloud we do not find direct correlations between Xmax and Ry. 

The positive slope in consistent with recent developments in theory and in- 
dicating alignment driven by the radiation field. The present data cannot con- 
clusively differentiate between direct radiative torques and alignment driven by 
H2 formation. However, the small values of A„iaa;(Ay=0), seen in several clouds, 
suggest a role for the latter, at least at the cloud surfaces. 
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The scatter in the Xmax vs. Ay relation is found to be associated with the 
characteristics of the embedded Young Stellar Objects (YSO) in the clouds. We 
propose that this is partially due to locally increased plasma damping of the 
grain rotation caused by X-rays from the YSOs. 

Subject headings: polarization, ISM: dust, ISM: Individual: Southern Coalsack, 
Chamaeleon I, Musca, p Ophiuchus, R CrA, Taurus. 



Introduction 



Starlight seen through the interstellar medium usually ends up being slightly polarized 
[at the level of up to a few percent), even when the background source is no t. This interstellar 
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interacting with the magnetic field. Theoretical models at tempting to explain inters tellar 



polarization soon followed with the first quantitative one by iDavis fc Greensteiru (1195 ll ) who 
identified grain alignment by the magnetic field through paramagnetic dissipation, as the 
origin of the polarization. This identification of dichroic extinction as the origin of the 
polarization has remained as the prime candidate up until this time, although the details 
of t he mechanism h as been modified by many authors over t he intervening half century 



see 



LazarianI (120031 ) for a recent review). For instance, while iDavis fc Greensteiru (Il95ll ) 
assumed paramagnetic grains, spun up by gas-grain collisions, it is now recognized that this 
combination will not suffice to keep the grains spun up against the damping influence of 
those same gas-grain collisio ns. The grain s have to have either a stronger magnetic moment 
(super-paramagnetic grains jMathislll986l )) or be spun up to rotatio nal velocities well above 
the thermal energy of the gas (suprathermal rotation (jPurcelllll979l )). 



The magnetic relaxation causing the rapidly spinning grain to align its axis of angular 
momentum with the local directi on of the magnetic fleld is well underst ood and based on 
solid state and nuclear physics (jPurcelll Il979l : iLazarian fc Draind Il999l ). The remaining 
poorly understood link in the explanatory chain to understand the origin of interstellar 
polarization is the mechanism of grain spin-up. Several mechanisms have been proposed on 
theoretical grounds, includin g the energy released from molecular hydrogen formation on 



the s urface of the dust grain (Puree 



fleld (jPolginov &: Mitrofanov 



1976 



19791) and direct torques from an anisotropic radi ation 



Weingartner k. Draind l2003l : ICho k. LazarianI l2005l ). 



It is generally agreed that, for the large grains involved in the polarization, the dampiri g 
of the grain spin, in neural gas, is due to gas-grain collisions ( iDraine fc LazarianI Il998l ). 
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For smaller grains or in regions with more extreme conditions , radi ation from the grains 
or other collisional partners can dominate (jPraine fc LazarianI Il998l ). Even so, the grain 
neutral damping is often used as a benchmark. A damping time can be defined such that 
it corresponds to the time it takes for a grain to collide with its own mass in gas particles. 
If we define the effective rad ius, a, as that which a grain of a given volume would have if it 
were spherical we find that (IWhittetlll992l ): 



td, 



lamping 



OC 



nyT 



T. 



where n and T are the gas number density and temperature, respectively. Hence, smaller 
grains have their rotation damped out more rapidly, and are therefore the hardest to keep 
spinning. A hotter and/or denser gas also dampens the grain rotation faster. Therefore, the 
size distribution of aligned grains can be used as a probe of the alignment mechanism. 

The first att e mpts a t determining the w'avelength dependen c e of i nterstellar polariza- 
tion (e.g. iHiltnen (jl949al )) were inconclusive. iDavis fc GreensteinI ( 1l95ll ) made a qualitative 
predic t ion o f the wavelength dependence of interstellar polarization, later expanded on by 
David (I1959I ). Isolated single wavelength observations hinting at such variations were indi- 
rectly reported by Stromgren in the 1954-1955 annual report of the Yerkes and McDonald 
Observatories where he notes that Hiltner had made observations at 1 and 2 /im, showing a 
steep drop in polarization compared to "the blue" (IStromgreru Il956l ) . The first systematic 
attempt at measuring the "polarization curve" to be published w as undertaken by A. Behr 
at Gottingen in 1958, at the prodding of L. Davis, Jr. (jBehrlll959l ). While wavelength depen- 
dence was indeed detected for a few stars, these observations only covered the wavelength 
range 0.37 - 0.51 fi m. The first obs ervations to span and locate the peak of the polarization 
curve were made by iGehreld ( 1l960l ) . 



When observed over the range of the optical and near-infrared spectrum, interstellar 
polarization takes on a characteristic wavelength dependence, which can be parameterized 
through the relation: 



exp[-K ln^{Xmax/>^)] 



(2) 



usually referred to as the Serk owski relation ( Serkowski 19731 ). if K is set to the fixed 
value 1.15, or the Wilking relation (IWilking et al.lll980l ) if K is used as a fitting parameter. 
Codina-Landaberry fc Magalhaes Jl976h first sugg e sted t hat K should b e rela ted to the 



size 



of the dust grains and results from Iwilking et al.l Jl982h and Iwhittetl Jl99jL 

indicate that 

K and Xmax are likely correlated and thus not independent parameters. (See IWhittetl (119921 ) 
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for an excellent review of the development of this parameterization) 



As has been shown by iKim &: Martini (119941 . Il995l ) the shape and variability of the 
polarization curve can be understood in terms of the size distribution of aligned grains 
and thus Xmnx provides a measure of the average size of the aligned grains. As noted by 



Kim fc Martini (119951 ) only a very small fraction of the grains are likely to be aligned and 
hence a degeneracy exists when using Xmax as a probe of grain alignment between the total 
grain size distribution and the fractional alignment in each size bin. This degeneracy is 
evident from the observational correlation of X max with Rv, the ratio of total-to-selective 



extinc tion, first noted by ISerkowskil (119681 ) (cf ISerkowski et al.l (Il975r): 



Whittet &: van Breda 



(119781 )). and is to be expected since, as shown by iKim et al.l (119941 ). R y traces the size 
distribution of the total g rain populatio n . As has been shown by e.g. IVrba fc Rydgren 



( 119841 ): IVrba et al.l (119931 ): IWhittet et al.l (120011 ). Ry is in turn, in general, correlated with 
the visual extinction, Ay. This is usually interpreted as being due to gra in coagulation 



causing the average grain size to increase at large r depths into the cloud (IBernard et al. 
I993I : Iwhittet et allboOll : IWurm k Schnaitei]l2002h . 



To break these degeneracies it is important to measure each parameter to high precision 
and to seek out regions where the grain size distribution and alignment conditions vary 
over relatively small scales. The latter should, particularly if radiative processes drive the 
alignment, most likely be found in the outer parts of the clouds where the radiation field seen 
by the grains varies relatively rapidly. For the former, we need regions where sightlines at 
similar observed visual extinctions show significant differences in the grain size distributions. 

A further complication in studying the variation of polarization is introduced by possible 
changes in strength or orientation of the magnetic field. To mitigate this concern it is 
important to use line of sight samples reliably constrained to probe single cloud. 

While the visual extinction is a convenient and straightforward probe of the amount of 
material along the line of sight, it is important to remember that the line-of-sight extinction 
may not be a good probe of the relative radiation field seen by the material being probed, 
particularly at modest values of Ay. The three dimensional geometry of the interstellar cloud, 
its relation to surrounding clouds and stars, and possible dumpiness of the material can cause 
the line-of-sight extinction to either over- or under-estimate the "effective extinction", here 
defined as the minimum opacity, vis-a-vis the diffuse interstellar radiation field, experienced 
by the material on the line of sight. As we will see below, this is an important consideration 
in interpreting optical polarimetry. 

The Southern Coalsack is a good target for studying the wavelength dependence of 
interstellar polarization and other absorption-based probes. The cloud is isolated and is well 
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located in three dimensions. Its location, straddling the Galactic mid-plane also guarantees 
a l arge number of backgro und sources, whether hot stars are required or not. As shown 



by lAndersson et al.l (120041 ) the outer parts of the cloud show a wide range of values of 
Ry. We have argued that this likely reflects clumpy dust destruction in the cloud envelope 
caused by the hot gas in the Upper Centaurus-Lupus Super bubble which envelopes the 
cloud. However, as we will discuss below, for emission tracers, and in particular continuum 
emission, the location of the Coalsack in the Galactic plane can be a major disadvantage, 
since foreground and background emission may be difficult, or impossible, to tell apart. 

For the purpose of the current study the important diffuse emission is the far-infrared 
light attributable to dust grains heated by the interstellar radiation field. Several authors 
have used color temperatures based on the IRAS 60/im-to-lOO/im ratio to show that observa- 
tionally this r atio p rovide s a tracer of the ra diation field impinging on the dus t clouds (e.g. 



Langer et all ([19891) (B5), ISnell et al.l (119891 ) (B18 and Heiles' Cloud 2) and iJarrett et al. 
Jl989h Ophiuchus)). While it is likely that the analyses in these studies are only qualita- 



tively valid due to the admixture of "large" and very small grains, and might not provide 
accurate absolute temperatures, the usage of the 60/xm-to-lOO/im ratio to trace the damping 



of the ISRF at different depths of a cloud is not in serious doubt (see e.g. iDraine fc Li 



(120071 )). In appendix |X] we use data from the clouds under study here to further support 



this usage. 

For those areas and line of sight where the interstellar extinction is dominated by a single 
cloud, or cloud complex (such as on high Galactic latitude sightlines) we can be relatively 
sure that the visual extinction and FIR emission are both caused by the same material. 
However, for low Galactic latitude clouds, where background emission is significant, and 
clouds with significant (high mass) star formation, where the internally generated radiation 
field is comparable to the ISRF, the correlation of the visual extinction and FIR emission 
can be expected to break down. 

In this study we have therefore complemented new multi-band polarimetry of the Coal- 
sack with archival polarimetry and other supporting data for five additional near-by inter- 
stellar clouds: Chamaeleon I (henceforth: Chamaeleon), Musca, p Ophiuchus (henceforth: 
Ophiuchus), R Corona Australis (henceforth: R CrA) and Taurus to address the alignment 
mechanism of interstellar grains. 

The remainder of this paper is organized as follows: First we present our new obser- 
vations of the Southern Coalsack. In section 13.11 we present and discuss the analysis of the 
polarimetry for both the new Coalsack data and the reanalysis of the polarization data for 
the other five clouds, and supporting data consisting of optical, near and far-infrared (FIR) 
photometry. We show that the wavelength of maximum polarization is correlated with the 
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visual extinction, albeit with many outliers. Section 13.21 proposes that the outliers in both 
the FIR and polarimetry data are due to cloud geometry and/or the presence of embedded 
- or nearby - stars making the observed visual extinction a poor tracer of the radiation field 
seen by the dust. Section [X^ shows that the FIR emission is anti-correlated with the visual 
extinction and that most of the outliers in the polarization plots can indeed be identified 
as outliers also in the FIR plots. This is expanded on in Appendix \M We then use the 
I(60/im)/I(100/im) vs. Ay relations to identify those sightlines where this is the case. Sec- 
tion 14.21 reanalyzes the Xmax vs. Ay relations for the different clouds after the screening 
performed in section 13.31 and we find that there are now tight correlations for four of the six 
clouds. No such correlations are seen for Xmax vs. Ry. However, the y-axis intercept for 
Xmax at Ay=0 is Correlated with < Ry >■ Section H73l discusses the dispersion seen in the 
Xmax vs. Ay relation in terms of the characteristics of the embedded young stellar objects 
in the clouds. 



Observations and Data Reduction 



We used the University of Cape Town Polarimeter (UCTP: ICropperl (119851 )) on the 1.9m 
telescopes of the South African Astronomical Observatory (SAAO), during 2005 April 6-8 
& 13 to perform multi-band observations of interstellar polarization of stars background to 
the Southern Coalsack. The UCTP was configured in the simultaneous linear and circular 
polarimetry mode, with a RCA31034A GaAs photomultiplier as the detector Measurements 
were perfo rmed using UBVRI fi lters in the Kron-Cousin system. We used HD 94851 and 
HP 98161 JTurnshek et al1 ll99oh as unpolarized stand ard stars and HD 155197, HD 154445 
( ISchmidt et al.lll992l ). HD 93632 (IMarraco et al.lll993l ) as polarized standard stars. At least 
one star from each class was observed each night. Measurements of the sky polarizations were 
acquired for each star immediately prior to an d followin g the main observation. The data 



were reduced using a custom software package ( Cropper! Il985l ). None of the stars showed a 
statistically significant circular polarization in any band. Target star information is given in 
Table [TJ Calibrated polarization for the Coalsack targets are summarized in Table [2] 



3. Analysis 

3.1. Polarization and Spectro-Photometric Data 



For each target in our sample we fitted the po larization dat a with both the Serkowski 
and Wilking relations. We then performed a F-test (lLuptonlll993l ) to determine whether the 
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additional parameter associated with the "Wilking" equation was justified. Only for those 
data sets where the additional parameter was justified at the 90% level did we accept it. For 
those stars where it was not statistically justified we have left the column for "K" blank in 
the result tables. The parameter uncertainties w ere calculated by t he fitting routines and 
verified using Monte Carlo simulations of the fits (jPress et al.l (119861 ). p529ff). Importantly 
for the present study, the wavelength of maximum polarization was rarely affected beyond 
the one-sigma level by the assumed or fitted value of K. This w as also the cas e if we instead 



of a fixed K-value of 1.15 used the relation K=1.66Amax+0.01 (jWhittetlll992l ). 



For all stars, the extinction parameters were recalculated and verified. Vis ual photom- 
etry was, where available, extracted from the Tycho database ( lH0g et al.ll2000l ) while Near 
Infrared (NIR) photometry was extracted from the 2MASS survey, except where explicitly 
noted. S pectra l classification was checked and updated. Intrinsic colors were estimated 



from lCoxl (120001 ) with uncertainties based on the reported uncertainties in the sources of the 
spectral classification and, where required, linear interpolation between the table entries. 



3.1.1. Coalsack Data 



Table [3] list the best-fit parameters for the polarization curve fits for the Coalsack stars. 
Since we do not have any Near Infrared polarimetry for the Coalsack stars, special care is 
needed to potentially justify the use of Wilking fits. For these stars we show the reduced 
for the fits. We note that while the F-test justifies the additional free parameter of the 
Wilking fit for several stars, in most cases the reduced fo^' these fits are then less than 
unity. The data and best-fit polarization curves are shown in Figure [TJ 

For HD 110432, the 2M ASS data is fiagged as being o f poor quality, we the r efore used 
an average o f the r esults from IWhittet &: van Bredal (Il980l ) , iDachs fc Wamstekerl (119821 ) and 
Dachs et al.l (119881 ) for this star. In Figure [2] we plot the location of the derived wavelength 
of maximum polarization (Xmax) as a function of visual extinction. As noted above, earlier 
studies of the dependence of Xmax on extinction parameters have found a rela tionship with 
the value of total-to-selective extinction (Ry; e.g. IWhittet fc van Bredal (119781 )). We do not 
find such a correlation for the Coalsack data. In Figure [21 we have color-coded the sightlines 
according to their Ry values. While the two sightlines above Ay=2.5 show Ry consistent of 
the value found at large visual extinction in this cloud (Ry ~3.25), no systematic trends in 
Ry are evident at smaller Ay. 



One possible caveat to the accuracy of the calculated Ry values comes from the fact 
that several of the target stars are in some spectral classifications listed as either emission 
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line stars or possible binaries. To evaluate the possible influence of these complications 
for the Coalsack sample we used the J-H vs. H-K diagram, plotted in Figure [31 Since 
the d irection of the red dening vector in this color-color diagram varies between regions 
(e.g. iRacca et al.l (120021 )) we calculated the H-K color excess that, for each target gave 



the smallest offsets between measured colors and best-flt reddened colors, given the spectral 
classiflcations of the stars and direction of the r eddening vector. We used ^jj^h /^h-k ~1.57 
( IKenvon et al.l 119981 ) and ^j_h/^h-k ~1.91 (INaoi et al.l l2006l ) as test cases. In Figure HI 
we plot the H-K color excess versus the derived visual extinction. The measured difference 
(H-K)-(H-K)o is plotted a x:es with error bars. Open symbols represent the color excesses 
calculated assuming 'Ej_h/'^h-k ~1.57, while flUed diamonds represent the color excesses 
calculated assuming 'Ej^h/'^h-k ~1.91. Only for HD 112661 and HD 112045 are offsets 
of more than 2a seen between the measured and best-flt H-K color excesses (2.1 and 2.2cr, 
respectively for 'Ejj_h/'Eh-k ~1.91). While both of these stars are listed in SIMBAD as 
multiple, neither shows an exceptional Yiy value. Even if these sightlines are suppressed, no 
correlation between Xmax and Rv is seen. 



3.1.2. Chamaeleon, Musca, Ophiuchus, R CrA and Taurus Cloud Data 



Figure [2] shows a tight correlation of X^ax with Ay over the limited range of Ay ~1.0-2.5 
but with outliers both a t sma ller and larger values of Ay- A similar weak correlation was 
noted by IWhittet et al.l (120011 ) in their Taurus data, but with a more pronounced scatter. 
To investigate whether the outlier points in the Coalsack plot are truly outliers in a real 
correlation, or whether the perceived correlation is instead a statistical fluke, we searched 
the literature for high quality, multi-wavelength polarimetry in sightlines penetrating well 
deflned interstellar clouds. 

We selected flve additional well-studied clouds with high quality published multi-wavelength 
polari met ry. Polarimet r y dat a were ex tracted from the st udies by: Ch amaeleon ; IWhittet 



(I1992I) fc ICovino et al.l (Il997l). Musca; lArnal et al.l (I1993I), Ophiuchus: Iwhitt^ Jl99i) fc 
Vrba et all (Il993h . R CrA: Iwhittetl dim and Taurus: Iwhittetl (Il992h . lwhittet et alTEoOli ). 



To ensure as uniform a dataset as possible in term of photometry and stellar characteristics, 
we extracted visual (B,V) and near-infrared data from the Tycho and 2MASS archives. We 
verifled, and where possible updated the spectral classiflcation, and then assigned intrinsic 
colors as for the Coalsack stars. The resultant stellar characteristics are listed in Table HI 
In a couple of cases unphysical values of - in particular - Ry are encountered, likely indi- 
cating problems with the spectral classiflcation. These sightlines were therefore excluded 
from the subsequent analysis. In Figure [51 we compare the values of Ry and Ay derived here 
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with those extracted from the hter ature. For both pa rameters, good agreement is seen. We 



note that the Ry values quoted in lArnal et aLl (119931 ) are based on polarimetry and not on 



photometry and we have therefore not included them in these plots. 

As for the Coalsack sightlines, we fitted Serkowski or Wilking functions to the polar- 
ization data, but selected only those stars with polarization measurements in at least 4 



f ilters (for this reason we also chose not to include the Chamaeleon data from lWhittet et al. 



(I1994J )). The polarization fit parameters are given in Table O In most cases the Xmax de- 
rived here agree with those in the original papers within la of the mutual uncertainties and 
in all but one cases within 2a. The one exce ption is the sight line towards HD 107875 in 



Musca (star number 7 in the nomenclature of lArnal et al.l (Il993l )). where the original paper 



reports A,„aa=0.641±0.009 /xm while we find Amaa;=0.49±0.01 /xm using their stated polar- 
ization curve parametrization. Given this discrepancy we have excluded this sightline form 
the analysis. 

Plotting Xmax vs. Ay foY all six clouds (Figure[6]), we see that, particularly for Chamaeleon, 
Musca and Taurus, a very similar structure is evident as for the Coalsack with a main group- 
ing of point lying along what seems like a linear correlation from Ay ~1 to ~2.5 mag. but 
with outliers both at small and large Ay. A critical issue, then, is whether the points we have 
here designated "outhers" are indeed that in a statistical sense, or whether, the suggested 
correlation of Xmax with Ay over approximately 1 to 3 magnitudes of extinction is illusory. 

As noted above, two independent ways exist for estimating the extinction for a given 
parcel of gas and dust. However, neither the directly measured visual extinction, nor the 
color-temperature of the FIR emission from the heated grains are immune to biases. We will 
next discuss how the combination of the two tracers of extinction can be used to mitigate 
biasing in the determination of the effective opacity. 



3.2. Effective extinction vs. line-of-sight extinction 

The line-of-sight extinction might not be a good indicator of the extinction - or equiv- 
alently, radiation field - seen by the material sampled by an absorption (line or continuum) 
experiment - the "effective extinction". Figure [7] illustrates some of the ways in which such 
discrepancies might occur. In this cartoon, we show a prolate cloud pointing towards the 
observer next to a roughly spherical cloud, with a "bridge" of material linking the two. The 
gray zone is meant to illustrate the part of the cloud into which the interstellar radiation 
field penetrates. 

For sightline 'A', which passes through the outskirts of the prolate cloud, the line of 
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sight extinction is much larger than the effective extinction seen by the material probed. The 
average radiation field seen by different parcels along the chord is similar and hence we'd 
expect this kind of sightline to show a large Ay while retaining a relatively large 60/100 
ratio. A similar effect might arise if the sightline passes within the sphere of influence of an 
embedded, or nearby, stellar source contributing to the radiation field (A'). 

For sightline 'B' the line of sight extinction, again, is larger than that seen by the 
material, but here, the different parcels along the chord see very different radiation fields. 
For this kind of sightline we would again expect a large Ay but here - since the FIR radiation 
traces booth radiatively heated dust and dust in the dark part of the cloud - expect a low 
60/100 ratio. 

For sightline 'C, which passes through a region between clouds (or in an inter-clump 
region of a clumpy medium) the line of sight extinction instead underestimates the effective 
extinction experienced by the material due to shadowing effects by the surrounding clouds. 

Finally, in sightline 'D' the observed and effective extinctions are similar and the different 
parcels of material along the chord experience similar radiation fields. 

Similar effects and observational discrepancies might also arise in a clumpy medium 
with the sightlines passing through predominantly clump or inter-clump material. 



3.3. Far Infrared Data 

Empirically an anti-correlation is seen between column density and dust t emperature 



based on the 60-to-lOO fim ratio in interstellar clouds (see, for instance figure 8 of ISnell et al. 



(119891 )). Since this ratio therefore indirectly traces the radiation field seen by the dust, we 
can use it as a probe of the "effective extinction" seen by the grain. While this emission 
is also likely prone to biases, such as variations in the fraction of small grains, we show in 
Appendix |X] that we can use the combination of visual extinction and Far Infrared (FIR) 
emission to screen the data sets for anomalous sightlines. We note that in this study we 
only do this in a relative sense for each cloud. We are not attempting to find absolute dust 
temperatures. While we will limit our analysis to the 60-to-lOO /zm ratio, we will use the 
25-to-lOO fim and 12-to-lOO ratios to argue in Appendix Rl that the local changes seen 
in the relation between the 60-to-lOO fim ratio and Ay are indeed likely due to irradiation 
differences. 



We extracted FIR data from the IRIS re-reduction (IMiville-Deschenes &: Lagachell2005l ) 
of the IRAS all-sky photometry. The spatial resolution of the IRIS maps are 3.8, 3.8, 4.0 & 
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4.2 arc-minutes respectively for the 12, 25, 60 and 100/im bands, but the maps are pixehzed 
on a 1.5' scale and hence over-sampled by about a factor of 2.5-2.8. For the current study 
we used a 3x3 pixel average that allows us to both lessen the impact of any emission from 
the background star and to identify it, if significant. Comparing to the IRAS point (faint) 
source catalog, we find that, in most cases, the stars - when at all detectable - make only a 
minor contribution to the total FIR light on our lines of sight. This is especially true of the 
longer wavelength bands. 

Since the nominal, single pixel, photometric u ncertainty in the IRIS maps i s mu ch 



smaller than the typical pixel-to-pixel variations (cf. iMiville-Deschenes fc Lagachd (120051 ) ). 
we found that the standard deviation in the 3x3 pixel averages can be used as a sensitive 
indicator of stellar contamination of the FIR flux densities, complementing the IRAS point 
source and faint source catalogs. For those sightlines where the background star - or a nearby 
star - contributes to the average, the standard deviation increases dramatically compared to 
the norm. 



3.3.1. Nominal Ay vs. I(60^m) /I (lOO^m) relations 

We can probe the difference between line of sight and effective extinction by comparing 
the observed visual extinction to the amount of emission due to heated grain, as probed by 
the ratio of 60/im to lOO/im flux densities. To find a nominal relationship for each cloud, we 
used sightlines with reliably determined visual extinctions and compared these to the FIR 
data. We used two partially overlapping data sets for this analysis: Field stars were first 
selected from the Hipparcos catalog, providing trigonometric parallaxes and thus allowing a 
clean separation of the stars into groups foreground and background to the cloud, but with a 
fairly limited total num ber of sightlines. W e also used the catalog of Tycho stars with known 



spectral classifications (IWright et al.ll2003l ) to maximize the number of sightlines used. In 
this case, we usually do not have explicit distance information and we therefore imposed 
somewhat more stringent selection criteria for which stars to include in the analysis. 

For the first field-star samples we selected all stars in the Hipparcos catalog within a 
three degree radius of /,6=(297,-15.5), (354, 15), (0,-19.5), (174,-14) and for the Chamaeleon, 
Ophiuchus, R CrA and Taurus clouds respectively. For the Musca cloud we use d a two degree 



radius centered on /,6=(301, -8). For the Coalsack, we used the target hst of ISeidensticker 



( 119891 ) as input in our Hipparcos search. From these original lists we then selected stars at 
distances (based on the Hipparcos trigonometric parallaxes) greater than those estimated 
for the clouds. Complementing the Tycho photometry with 2MASS NIR photometry we 
calculated visual extinctions. We rejected stars with negligible extinction and were left with 
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samples of 39, 30, 23, 34 and 27 field stars for Chamaeleon, Musca, Taurus, R CrA and 
Ophiuchus clouds respectively. 



For the second sample we selected stars as above from the catalog of IWright et al. 



( I2OO3I ) for Chamaeleon, Coalsack, Musca, Ophiuchus and R CrA. For Taurus - because of the 
elongated shape of the cloud - we selected stars within a three degree radius of two centers 
at i?^,L'ec= (04:40:00,25:30:00) and (04:15:00,28 00 00). As with the Hipparcos sample, 
we extracted Tycho and 2MASS photometry, and calculated visual extinctions and ratios 
of total-to-selective extinctions. Based on Ay vs. distance diagrams from the Hipparcos 
samples, we then selected only those stars with Ay >{0.3, 0.4, 0.4, 0.4, 0.3 and 0.5} for 
Chamaeleon, Coalsack, Musca, Ophiuchus, R CrA and Taurus respectively. We used the 
IRAS point- and faint-source catalogs to screen out stars detected as point sources and finally 
screened for unreasonable values of Ry, which likely reflect unreliable spectral classifications. 
Finally, we eliminated small regions (usually 0.5 or 1 degree radii) around stars where the 
Iqo/Iioo ratio showed locahzed influence from those stars (see Appendix |A]) . We were then 
left with 90, 231, 60, 58, 96 and 154 stars respectively for the six clouds. 

Because of the relatively bright limiting magnitude of the Hipparcos database, the size 
of the area used to select the field stars involves a trade-off between the largest accept- 
able distance from the cloud center and number statistics. While the polarization sam- 
ples do cover smaller areas on the sky, the three-degree ra dius (two degrees for Musca) 



was chosen based on emission-line and extinction maps (e.g. iCambresyl (119991 )). Although 
we do not have kinematic data tracing the material giving rise to the extinction in the 
field star samples, we can use the measured Ay values to see that this sample likely does 
probe the outer layers of the molecular clouds under study. Our Hipparcos field star sam- 
ples have visual extinctions of Ay{min, max,mean}={0.3,2.4,1.0} mag. for Chamaeleon, 
{0.1,2.4,0.9} mag., for the Coalsack, {0.4,1.7,0.77} mag. for Musca, {0.4,4.5,1.6} mag. for 
Ophiuchus, jo. 1,1. 7 ,0.6) mag, for R CrA and {0.4,3.3,1.3} mag. for Taurus. As shown by 
Savage et al.l (119771 ) and iRachford et al.l (120021 ) the transition to molecular hydrogen occurs 
already at E^-y ~0.1 (Ay ~0.3) and hen ce most of our field sta r sightlines probe mo l ecular 
mat erial. In addition, as shown by e.g. IWannier et al.l (119831 ) . Ivan der Werf et al.l (119891 ) 
and Andersson fc Wanniei (Il993), raolecul ar clouds are surrounded by extensive atomic en- 
velopes. Figure 3 of iBoulanger et al.l (119981 ) also illustrates this for the Chamaeleon complex. 
Particularly for the high-latitude clouds in our study, it is thus unlikely that the field stars 
probe material unrelated to the clouds probes by the polarimetry. 

As we are explicitly searching out, and expecting to find, points w here the sys t emati c 
errors dominate the random ones, we used a robust fitting algorithm (jPress et al.l (119861 ) . 
p 539ff) to find the nominal I(60/im)/I(100/xm) vs. Ay relations for each cloud. This 
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algorithm uses an iterative procedure based on Tukey's Biweight weighting with a hmit of 6 



outher-resistant standard deviations ( Press et a. 



19861 ). as implemented in the IDL routine 



"robustjinefit", available in the "astron" librarja. For the larger Tycho sample we also used 
weighted linear fits. The two algorithms yield very similar fitting parameters in this case. 

For the Hipparcos samples, the robust fits for Chamaeleon, Musca, Ophiuchus, R CrA 
and Taurus yield: 



(0.241 ± 0.005) (0.022 ± 0.005) [Cham] 

(0.29 ±0.01) (0.026 ± 0.010) [Mmsc] 

ho/hoo = (0.27 ± 0.02) - Ay x (0.01 ± 0.02) [Oph] (3) 

(0.21 ± 0.01) (0.035 ± 0.01) [R CrA] 

(0.196 ± 0.008) (0.013 ± 0.006) [Tau] 

For the Tycho samples, the robust fits yield: 

(0.240 ±0.004) (0.018 ±0.004) [Cham] 

(0.29 ±0.1) (0.030 ± 0.007) [Mmsc] 

ho/hoo= (0.27 ±0.1) - Ay x (0.024 ± 0.010) [Op/i] (4) 

(0.244 ± 0.006) (0.023 ± 0.007) [R CrA] 

(0.179 ±0.003) (0.012 ± 0.002) [Taw] 



In Figure [8] we show the best linear fits (solid lines, top panel) using the robust fitting 
algorithm to the Tycho field star samples only (i.e. the polarization targets were not included 
in the fits). 

For Chamaeleon, Musca, Ophiuchus, R CrA and Taurus well-defined linear correlations 
are found (for Ophiuchus, the uncertainties on the fit coefficients for the Hipparcos sample, 
are too large for the fit to be significant). For the Coalsack, which straddles the Galactic 
plane, this technique fails to provide a reasonable correlation, presumably because of the 
strong influence of diffuse background sources in the FIR data. 

The middle panels of Figure [8] shows the locations of the polarization targets in the FIR 
vs. Av diagrams. 

In the lower panels of Figure [HI we show the distributions of offsets in leo/Iioo from 
each data point to the best fit line at the same value of Ay (dashed histograms for the field 
stars and solid histograms for the polarization targets), and the best fit Gaussians to these 



^Available at |http: / /idlastro.gsfc.nasa.gov7| 
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distributions. The widths of the distributions are similar for all clouds (cr=0.012, 0.014, 
0.036, 0.026 and 0.016, respectively). Based on similar geometrical arguments as for Figure 
[3 it is likely that these widths trace the porosity of the outer parts of the clouds where the 
radiative grain heating occurs. In Appendix |X] we show that the widths of these dispersions 
are correlated with the SFR in the clouds. 

The dashed lines in Figure [8] show the 2a distance from the best-fit solution. Points 
below and to the left of the left-hand dashed lines are likely to have characteristics of sightlines 
'B' in Figure [71 for high values of Ay, or 'C, for low values of Ay, while points above and 
to the right of the right-hand dashed lines are likely to have characteristics of sightlines 'A'. 

We note that with the exception of the identified outliers, the polarization sample targets 
and the field stars overlap in the well-defined Gaussian distributions in offset from the fits, 
providing additional support for the assumption that the two observational samples are 
drawn from a common parent population. 

While the above screening method fails for the Coalsack, inspecting the \max vs. Ay 
plot for the Coalsack and comparing it to those of Chamaeleon and Taurus, we can see 
that the points at relatively high extinction and relatively low Xmax in the Coalsack have 
corresponding sightlines in both Chamaeleon and Taurus and that these sightlines in the 
Coalsack fall in the equivalent "A" region in the leo/Iioo vs. Ay plot as their counterparts 
for Chamaeleon and Taurus. Similarly, the points at high Xmax and low Ay are similar to 
points in region "C" for Chamaeleon both in the \max vs. Ay and leo/Iioo vs. Ay plots. 
While based on a more indirect and less satisfactory procedure, we therefore screen these 
points also in the Coalsack sample. In Figure [9] we show the resulting leo/Iioo vs. Ay plot 
with the sources thus deselected marked. The width of the Aleo/Iioo distribution for the 
remaining sightlines for the Coalsack is (T=0.04 also similar to the values found for the other 
clouds. 



4. Results and Discussion 

4.1. Origin of the Outliers 

We argue that the outliers in figure [8] are likely primarily due to localized differences 
in the irradiation of the dust. Appendix lA] provides several direct lines of evidence in favor 
of this interpretation. For Chamaeleon and R CrA all the sightlines identified as type "A" 
can be seen to be due to proximity to HD 97300 or IID175362, respectively. Similarly, for 
Ophiuchus almost all type "A" sightlines can be seen to be close to either o Sco or p Oph D 
(Figure [To]). Additionally, as also shown in appendix [Xj the ratio leo/Iioo shows a stronger 
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response to the vicinity of a hot star than I25/I100 while the I12/I100 ratio shown httle or no 
response to the proximity of hot stars. These results provide strong support for an origin in 
irradiation differences. As shown in Appendix |X] the dispersion around the best-fit lines in 
the Igo/Iioo vs. Ay plots are correlated with the characteristics of the star formation in the 
cloud (see below). We interpret this as due to increasing porosity (dumpiness) in the clouds 
produced by increasing star formation activity, which changes the radiative transfer of the 
light heating the dust grains. A quantitative analysis of this clumpy transfer is beyond the 
scope of the present paper. 

It should be noted that the screening procedure used in figure [8] does not rely on 
this interpretation. The screening as such is a straightforward numerical procedure and 
uses a explicit numerical limit for which sightlines to label "anomalous". The use of the 
I(60/im)/I(100/i) vs. Ay ratio to screen the polarimetry data should also not introduce any 
biases in the latter. This is as on ly a very small fraction of the grains contribute to the 
polarization (IKim fc Martini 119951 ) while all grains contribute to the FIR emission. Also, 
even if the outliers in figure [8] were due to significant enhancements in the population of 
very small grains, the size distributions of the VSGs and the grains responsible for the visual 
polarization have, at most, a minute overlap (IKim fc Martini 1 1 9951 : iDesert et al.lll99Cll ). 



4.2. Debiased Xmax vs. Ay plots 

In Figure [11] we show the plots of Xmax vs. Ay which result if we reject the points in 
Figure [8] which falls more than 2a from the robust fits for the Tycho sample. Overlaid are 
the best linear fits for each cloud (full drawn lines). The dashed lines represent the relations 
based on the average slope from all the clouds and a zero intercept based on the mean value 
of Ry of the cloud (see below). Very similar plots result for the Hipparcos based screening, 
with the main difference that the R CrA sample shows a broader scatter. We show the Xmax 
vs. Ay plot for R CrA using the Hipparcos screening in Figure [131 In what follows, we will, 
except were explicitly noted, discuss these "de-biased" data sets. 

Three conclusions can immediately be drawn from these new plots: 

1) While the plot of Xmax vs. Ay including all the data shows only a very weak cor- 
relation, the "screened" plot shows a distinct correlation of Xmax with Ay for Chamaeleon, 
Coalsack, Musca and Taurus and consistent results for Ophiuchus and R CrA. We can quan- 
tify this latter statement by performing Spearman rank order tests on the two groups of 
data sets. For Chamaeleon and Taurus both the unscreened and Hipparcos screened data 
sets show small (<6%) probabilities for being uncorrelated (the Tycho screened sets yield 
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2% and 14% probabilities). For the Coalsack the probabihty for an accidental correlation 
drops from 44% to 6% after screening. For Musca (where no sightlines are screened) we find 
a 18% probability of being uncorrelated. 

2) Within each individual cloud sample, no obvious correlations are seen between Xmax 
and Ry. 

3) It now is clear to the eye that the slopes of Xmax vs. Ay in the individual clouds are 
very similar, but the three groups are offset from each other in Xmax- We used robust fits to 
quantify this and find that for the Hipparcos screened sample: 



(0.54 ±0.02) 


(0.04 ±0.01) 


[Cham] 


(0.50 ±0.02) 


(0.05 ±0.01) 


[CS] 


(0.55 ±0.02) 

(0.60 ±0.05) ^ 


(0.04 ±0.02) 


[Muse] 


(0.01 ±0.03) 


[Oph] 


(0.72 ±0.10) 


(0.01 ±0.08) 


[R CrA] 


(0.48 ±0.02) 


(0.05 ±0.01) 


[Tau] 



(5) 



While for the Tycho screened sample: 



(0.54 ±0.01) 
(0.50 ±0.02) 
(0.55 ±0.02) 
(0.64 ±0.05) 
(0.75 ±0.03) 
(0.51 ±0.02) 



+ Avx 



(0.024 ±0.006) [Cham] 
(0.05 ±0.01) [CS] 
(0.04 ± 0.02) [Muse] 
(0.01 ± 0.03) [Oph] 
(0.03 ± 0.02) [R CrA] 
(0.03 ± 0.02) [Tau] 



(6) 



The slopes in Xmax vs. Ay are all very close and, in particular, those for Chamaeleon, 
Coalsack, Musca and Taurus are all both very close and have small error bars, indicating 
that the slope is universal (Figure [TH left panel). A weighted average of the slopes yields 
dA^ax/dAy=0.028±0.005 (0.038±0.007) for the Tycho (Hipparcos) screened samples. 



As has been shown by IWhittet fc van Bredal (119781 ) Xmax is generally found to be cor- 
related with Ry. Thi s is not surprising as, as has been shown by lKim &: Martini (119951 ) and 
Kim fc Martini (119941 ). changes in both Xmax and Ry are most sensitive to changes in the 



population (total and aligned part) of the smaller grains. In the right-hand panels of Figure 
[11] we plot Xmax vs. Ry for the screened samples. No correlations are seen even for the clouds 
where correlation in the Xmax vs. Ay plots are evident. 



We calculated weighted averages of Ry for the full polarization targets samples and find 
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Ry=3.6±0.1, 3.6±0.2, 3.6±0.2, 3.7±0.7, 4.1±0.2 and 3.3±0.1 for Chamaeleon, the Coalsack, 
Musca, R CrA and Taurus, respectively. In Figure [H] (right panel) we plot the result with 
the best linear fit of the zero-intercept of Xmax {Xmax{-^v=0)) for the Tycho screened samples 
with the average Ry. The best, weighted, fit yields: 

>^max{Av = 0) = (-0.47 ± 0.16) + (0.28 ± 0.04) < Ry > (7) 
or, if we impose a zero intercept: 

Ar„ax(^y =0) = (0.151 ± 0.002) <i?y> (8) 

for the Hipparcos screened sample the equivalent fits yield: 

A™ax(^y = 0) = (-0.26 ± 0.24) + (0.22 ± 0.07) < Ry > (9) 
or, if we impose a zero intercept: 

>^maxiAv = 0) = (0.146 ± 0.002) <Rv> (10) 
In both cases th e zero intercept fit is close to the Amaa;=(0.18±0.01)Rv found by 



Whittet fc van Bredal (Il978f ). 



We interpret the correlation of Xmaxi^v = 0) with < Ry > as being due to the under- 
lying differences in the total grain size distribution between the clouds, and the (universal) 
slope in Xmax vs. Ay as due to varying degrees of grain alignment at different depths into 
the clouds caused by a common mechanism. 

The data indicate that the lack of aligned small grains at larger visual extinctions, indi- 
cated by large values of Xmax, is due to loss of alignment rather than grain destructions. This 
can be seen by comparing Xmax as a function of Pmax/Ay and of Ry. In Figure [TS] we show 
Xmax for the Tycho screened lines of sight, adjusted to Chamaeleon, as functions of either 
the alignment efficiency (pmax/Av) or R^. The adjustment performed here consists of sub- 
tracting the difference in derived Xmaxi^v = 0) between the cloud and that for Chamaeleon 
for the sightlines in each cloud (i.e.: A;,,,=A^a,-(A^,,(Ay=0)'=''"^'^-A„,,(Ay=0)^'^'^"^.)). 

A distinct anti-correlation (albeit, again, with outliers) is seen in the Xmax vs. Pmax/Ay 
(Figure [T5| left panel), indicating that as the relative number of aligned grains increase the 
average size of the ahgned grains decrease. No correlation is seen in the Xmax vs. Ry plot 
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(Figure [T5|, right panel) Very similar plots, again, result if we instead use the Hipparcos 
screened sample. 

As equation [T] shows, when all other parameters remain fixed, the smallest grains will 
have their rotation damped out by gas collisions the fastest. Hence, when we find that smaller 
and smaller grains remain aligned either the damping is lessened or the driving mechanism 
for the spin-up in being enhanced. The correlations of Xmax with Ay shows that as we get 
increasingly close to the cloud surface, smaller grains remair i aligned. As has been shown 
on large scales for low radiation intensity cloud envelopes by lAndersson fc Wannierl (1l993l ): 
Wannier et al.l (119991 ) and recently on sma ll scales in the higher radiation field case of the 
Horse Head nebula by lHabart et al.l (120051 ) . cloud envelopes show isobaric structures in the 
gas pressure and hence equation [T] becomes 



td, 



ampmg OC 



gas 



with Pgas a constant for a given cloud. 



Hence an increase in temperature (assuming an isobaric equation of state) should lead 
to a decrease in Xmax- We searched the literature for Copernicus and FUSE measurements of 
the J=1/J=0 excitation temperature for the six clouds in our study. Unfortunately, only for 
Ophiuchus and Chamaeleon are there presently multiple interstellar sightlines published and 
not all of those sightlines have b e en stu died in polarimetry. For O phiuchus, five si g htline s 
were studied by ISnow fc Jenkins! (Il980l ) with an additional two bv ICartledge et al.l (120041 . 
These yield a range in 1-0 temperatures of 46 to 90 K. For Chamaeleon, Grv et al. ( 2002h 



report 1-0 excitation temperatures for three stars between 60 and 66 K. iRachford et al. 



(120011 ) report T=63 K for the line of sight towards HD110432 behind the C oalsack. For the 
five O phiuchus sighthnes with H2 data, measured Xmax can be extracted from lSerkowski et al. 



(I1975I ) while one of the Chamaeleon stars (HD 96675 = F24) and HD 110432 are included 
in the present data samples. In Figure [16] we plot both the measured values of Xmax and 
the offset from the best-fit relations in Xmax vs. Ay for each cloud. No clear correlation is 
apparent. We therefore conclude that variations in gas temperature are unlikely to explain 
the variations in the grain alignment. 

Based on these observational results we conclude that the data support grain alignment 
driven by the radiation field. This conclusion is consistent with recent developments in 



the theory of interstellar g r ain alignment. Sever al authors (jPraine fc Weingartnerl Il996 



Weingartner fc Draindl2003l : ICho fc Lazarianll2005l ) have shown that direct radiative torques 
are the theoretically most likely mechanism for explaining the spin-up of interstellar grains 
required to allow the magnetic alignment to take place. 
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We note, however, that our results cannot exclude alignment driven by molecular hy- 
drogen formation (in steady state any formation of H2 must be preceded by the radiative 
destruction of the molecule, and henc e this mechanism is als o more strongly active at smaller 
Ays). Indeed, if we use the models of iKim &: Martini (Il995l ) to estimate the smallest aligned 
grains at the surface of the clouds, using Amaa:(Av'=0)~ 0.5 fim, we find that grains as small 
as 0.01-0.04 fim need to be at least partially a li gned. For silicate grains in a diffuse cloud 
environment. Table 5 of iDraine fc Weingartnerl (119961 ) shows that for such small grains, H2 
for mation driving is dominan t over direct radiative torques. We note though that according 



to iLazarian fc Draind (119991 ) very small grains should flip frequently due to Barnett and 
nuclear relaxation. Thus under H2 formation torques (fixed in the grain coordinate system) 
they are expected to become thermally trapped and not be able to achieve suprathermal 
spins. Further observational and theoretical studies are clearly needed to clarify the origin 
of the very small Amax(Ay=0) we observe. 



4.3. Influence of star formation 

Figure [TT] indicates that the scatter in Xmax is related to the star formation activity of the 
cloud. We will here explore possible mechanism for this dependence, but start by noting that 
turbulence, and its associated line of sight variations in the magnetic field direction, is not a 
likely cause. This is because Xmax is not, to first order, dependent on the absolute amount of 
polarization and since turbulence should affect the polarization in all bands similarly, leaving 
Xmax unaffected. 

In Figure [T7| we show histograms of the distance from the best fits in Xmax with Ay for 
the six clouds using the "universal slope" relations. As is clear by a visual inspection, the 
scatter increases from the Coalsack, Chamaeleon and Musca through Taurus to Ophiuchus 
and R CrA. 



In Figure [TS] we plot the scatter around the best fit in Xmax with Ay as a function of 
the "bright star fraction" (bsf(L^g;)) of embedded Young Stellar Objects (YSO). The bright 
star fraction is defined as the number fraction o f objects brighter than a given bolometric 
luminosity (L^^^) to the total number of objects (cf lChen et al.l (119971 )). The error-bars on the 
bsf:s reflect counting statistic s. The data for Chamaeleon, Ophiuchus, R CrA and Taurus are 
taken from lChen et al.l (119971 ). For the Coalsack, which does not show any star formation, we 
here use bsf=0. Musca has not been specifically studied for star formation act ivity although 



some T Tauri stars are detected in the general Chamaeleon-Musca complex (IMizuno et al. 



19981 ). To separate the Coalsack and Musca points we have for plotting purposes assigned 
a bsf of 0.17 mid- way between Chamaeleon and the Coalsack. Neither point was used the 
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fitting. 



We used power law and exponential fits for several choices of the break-point (L^^^; see 
above) in the bsf and find a best fit using an exponential for Li,oi=11-'q, yielding a correlation 
coefficient of R=0.88 (for the Hipparcos screened sample and 0.79 for the Tycho screened 
sample) the resulting fit is over-plotted in Figure [181 At both smaller and larger L^^^ the 
fit is worse. While not conclusive, this is consistent with an origin of the correlation due 
to the infiuence of the YSOs since the contribution, per unit luminosity range, to the total 
luminosity of the YSOs peaks around 1-2 Lq. This is illustrated in Figure [T^ where we show 
the relative contribution to the bolometric and X-ray luminosity per unit Ltoi for different 
sub-samples. The data for the bolometric lumi nosity are a combination of the data for 
Chamaeleon, Ophiuchus and Taurus, taken from I Chen et al.l (Il995l ). The individual cloud 
data give sini i lar re sults but with bigger scatter. The data for the X-ray luminosity are from 
Grosso et al.l (120001 ). We chose a variable binning aimed at providing good resolution while 
collecting enough targets per bin to achieve reasonable statistic per bin. The relative binned 
bolometric luminosity peaks at ~1 Lq depending marginally on the exact binning, while 
the X-ray luminosity peaks at ~l-2 L0 with significant contributions from individual high 
luminosity stars. (The 10-20 L^y^ bin however only contains two stars. We have also excluded 
target A21 from lCrosso et all fl200oh with 1^=11000 and log(Lx)=30.8 in the plot.) 



In Taurus (and R CrA if using the Hipparcos screening) the outliers in Figure [T7| are 
primarily located on the positive side of the plots. Depending on whether we use the Tycho 
or Hipparcos based screenings, the data for R CrA either shows a fairly wide dispersion or a 
significant offset in the center of the distribution. For Chamaeleon, the Coalsack, Musca and 
Ophiuchus, the average of the distances of points from fits are much less than their associated 
standard deviation. Also for Taurus the average offset is smaller that t he standard deviation, 
but here the distribution has a statistically significant positive skew (jPress et al.l ( 119861 ) . p. 
457) of 0.66±0.53 (for the Hipparcos screening; 0.57±0.53 for the Tycho screening). 

The scatter is likely due to porosity in the clouds, introduced by the effects of star 
formation, as seen in the dispersions between leo/Iioo vs. Ay plots (Appendix |A|). However, 
for variations in the radiation field intensity caused by clouds porosity alone we would expect 
that the scatter be symmetrical. Particularly for Taurus and R CrA we are thus led to 
consider possible mechanisms which will drive Xmax selectively to larger values, and hence 
also possible sources of additional, localized, grain rotation damping. 

Given the ubiquity of X-ray emission from YSOs, it is worth considering what effects 
X-rays from the embedded YSOs would have on the grain rotation. One important X-ray 
induced effect is the ejection of photoelectrons from the gr ains and both the s ubseq uent 
heating of the gas and charging of the grains. As shown by iDraine &: LazarianI (119981 ) the 
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dominant rotation-damping mechanism in molecular cloud environments for very small grain 
(smaller than those discussed herein) is plasma drag due to the interaction between the ions 
in the gas and the electric dipole mom ent of the grain. Since the dipole moment of the 
grain is proportional to the grain charge (jPraine &: Lazarianlll998l ). enhanced localized grain 
charging would also mean enhanced localized rotation damping, in addition to the enhance- 
ment from gas hea t ing. If we extrapolate the molecular cloud (MC) panel of Figure 4 of 
Draine &: LazarianI (119981 ) to the 0.05 /xm range, we find that in the nominal model, neutral- 
gas drag is the dominant damping mechanism. However, plasma drag is only less important 
by a factor of a few. Hence, an increase in the grain charge by an order of magnitude would 
invert this relative importance. The typical X-ray (hi/~0.5-5 keV) lur ninosity of embed- 



ded YSOs is 10^^- 10^^ ergs s ^ with flares reaching significantly higher (jOzawa et al.l 12005 
Grosso et al.lboooh . 



To derive a first order estimate of whether X-rays from YSOs can significa ntly effect the 



grain rotation, we use the Ophiuchus cloud observations. In the observations of iGrosso et al. 



( I2OOOI ) for the inner part of cloud, 66 targets were found within an approximate radius 
of 15', or at a distance of 140 pc, 0.6 pc, yielding a characteristic projected source-to- 
source distance of 0.075 pc. At this radius the X- rays from a single s ource contributes 
an energy density of about u=5xl0~^^ ergs cm~^. iDraine &: LazarianI (119981 ) assumed a 



energy density in the radiation field of the 1% of the nominal ISR F for their molecular cloud 
medium, or u~9xl0"^^ ergs cm~^ (jWeingartner fc Draindl200lh . Taking into account also 
the contributions of F UV emission from the YSOs, the high photoelectric yield at soft X-rays 
JPwek fc Smithlll996h and the con tributions from the diffuse X-ray emission, likely caused 
by stellar winds (lEzoe et al.ll2006l ). we feel that the influence of embedded X-ray sources 
is a viable cause of the localized rotation damping and warrants further study. However, 
observational studies and modeling, beyond the scope of the present paper, are required to 
quantify the importance of X-ray induced grain rotation damping. 



5. Conclusions 

We have used new multi-band polarimetry of the Southern Coalsack to study the align- 
ment of interstellar grains. For this study the Coalsack has the advantage of not showing 
a systematic increase in the total-to-selective extinction, Ry with increasing visual extinc- 
tion and to show no star formation. The former allows us to break the degeneracy between 
increasing values of the wavelength of maximum polarization (Xmax) due to the potentially 
lessened driving of the grain spin and due to grain growth. We find a tight correlation of 
Xmax with Ay ovcr the limited range of Ay=l-2.5, but with a significant amount of outliers. 
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Using polarization data from the literature, together with archival FIR emission data, for 
an additional five near-by clouds, we show that the outliers in the Xmax vs. Ay plot are 
primarily due to geometrical effects, causing the measured visual extinction to become a 
poor tracer of the extinction "seen" by the dust grains under study. 

Screening the observations based on the leo/Iioo vs. Ay relationship we find that for 
clouds without, or with only low star formation activity, Xmax is tightly correlated with Ay. 
We performed the screening based on two somewhat overlapping data sets, namely field 
stars background to the clouds from, either the H ipparcos catalog, o r from the catalog of 



Tycho targets with known spectral classifications (IWright et al.ll2003l ). Screening with the 
two field-star samples yield consistent results but slightly different numerical parameters. 
In both cases, the correlation between Xmax and Ay shows a universal slope and zero-Ay 
intercepts proportional to the averages of Ry in the clouds, such that 

Xmax = (0.15 ± 0.01) <Rv> +(0.038 ± 0.007)Ay (12) 
for the Hipparcos based screening and 

Xmax = (0.15 ± 0.01) <Rv> +(0.028 ± 0.005)Ay (13) 

for the Tycho based screening, where we have imposed a zero intercept in the X^^^{< Ry >) 
relation. 



Within each cloud, we do not find correlations between Xmar and R 



V 



We interpret the positive slope in Xmax vs. Ay as evidence for radiation induced grain 
spin-up. Our data cannot conclusively differentiate between direct radiative torques and H2 
formation torques. A possible role f or the latter is howev er indicated by the small value of 



Xmax at Ay=0, usiug the models by iKim fc Martini (119951 ) 



Compensating for the different Xmax{-^v=^) between clouds, we find that Xmax is anti- 
correlation between with the alignment efficiency, Pmaa/Ay, further supporting the conclu- 
sion that the systematic variation in Xmax as a function of Ay is indeed due to grain alignment 
rather than changes in the grain size distribution. 

The scatter around the Xmax vs. Ay correlation - as well as the leo/Iioo vs. Ay relation- 
ship - increases with the star formation rate in the cloud and is correlated with the number 
fraction of YSOs brighter than 1 Lboi. In particular for the R CrA and Taurus clouds the 
scatter shows asymmetries toward positive values, indicating possible evidence for localized 
enhancements in the grain rotation damping. The implied patchiness of the excess damping 
together with the correlation of scatter in Xmax with bright YSOs and the high X-ray output 
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of such embedded sources lead us to propose a connection to photoelectric grain charging 
and associated plasma damping of the grain rotation. 

Thus, to quantitatively understand interstellar polarization we must know the strength 
of the radiation field at the location of the grain, as well as the star formation environment 
of the region, likely including the X-ray flux seen by the grain. Further surveys of the 
polarization curve for clouds with known star formation rates (SFR) as well as significant 
numbers of well-characterized line of sight at large (Ay > 3 mag.) visual extinctions, both 
in the present sample of clouds and in new ones would allow the results presented in this 
paper to be tested and extended. We are currently pursuing several such studies. 

It is a pleasure to acknowledge many helpful discussions with Alex Lazarian, Joe Wein- 
gartner and Bruce Draine. Kip Kuntz and Steve Drake provided helpful pointers to the 
X-ray properties of YSOs. The referees provided several suggestions and challenges that 
help to significantly improve the paper. This research has made extensive use of the VizieR 
catalogue access tool at CDS, Strasbourg, France, as well as of the NASA/ IPAC Infrared 
Science Archive, operated by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space Administration. B-G A. 
gratefully acknowledges travel support from the Center for Astrophysical Sciences at the 
Johns Hopkins University. 



A. Does the 60-to-lOO/im ratio really trace the local radiation field intensity? 



As noted in the introduction, several authors (e.g. 



Laneer et al. 


( 


1989): 


Snell et al. 



(1l989l ): IJarrett et al.l (1l989l )) have shown that the leo/Iioo ratio is anti-correlated with col- 



umn d^n^it3^_of_CO and interpreted this as a temperature effect. However, as shown by 



e.g 



Bernard et al.l (Il993l ) the relative abundance of the very small grains (VSG) decreases 
towards the center of clouds and we therefore need to ascertain what part of the variations 
seen in the FIR ratios is dominated by grain heating rather than by abu ndance effects. First, 
we note that modern theories of the FIR emission (IDraine fc Lil 120071 ) supports the use of 
the leo/Iioo ratio as a irradiation tracer. Figure 12 of that paper shows a very small effect 
from the variation of the PAH fraction, while figure 13 shows that the ratio is very sensitive 
to the strength of the radiation field. 



Arguing in somewhat more detail, the models by lDesert et al.l (Il990l ) indicate that the 
100 /xm band emission is dominated by large grains with a minor contribution from the very 
small grains (VSG), the 60 fim band contains a significant contribution from both large grains 
and the VSG. The 25 fim band emission predominantly samples the very small grains (with 
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some emission from poly-aromatic hydrocarbons; PAH) and the 12 fim band predominantly 
samples the PAH:es with a contribution from the VSGs. Hence, if the variability in the 
FIR ratios is primarily driven by VSG abundance effects, the I25/I100 ratio should show 
bigger deviations whereas if the variability is primarily due to temperature variations, the 
leo/Iioo should be more affected. Also, if irradiation (and hence heating) dominates the FIR 
emission, we should be able to detect effects due to physically near-by hot stars on the dust 
in each cloud. 

In Figure [20] we plot the ratios of leo/Iioo, I25/I100 and I12/I100 as functions of the on- 
the-sky distance between the locations of the Tycho selected field stars and early B stars 
located close to the clouds (Table [6]). This figure shows both that the leo/Iioo responds to 
the enhanced radiation field caused by localized sources and that the I25/I100 and I12/I100 
ratios show successively less effect from the enhanced radiation field. 

Figure [22] shows the dispersions around the best fit in leo /Iioo vs. Ay for the Tycho 
selected field stars for the clouds in our study as a function of the bright star fraction for each 
cloud (we have here also added the result from an equivalent analysis or the Lupus I field). 
A clear correlation is seen, tightly fit by a power law function. This can be qualitatively 
understood in terms of variations in the local radiative heating of the dust in a porous 
medium where the porosity is due to the influence of newly formed stars. 

Finally, we must ask whether the field star samples provide accurate representations 
of the behavior of the dust sampled by the polarization targets. The field star samples are 
generally somewhat shallower than the polarization samples. However, as the lower panels of 
Figure [8] shows - once anomalous sightlines are excluded - the distributions of the field stars 
and polarization targets around the best fit lines are very similar. Unfortunately, background 
stars with well established spectral classifications and large visual extinctions are very rare. 
However, in Figure [2T] we s how the Im/Imn vs. A y pl ot for Taurus, with t he fie ld stars 



from the water ice surveys of [Murakawa et al.[ (|2000[ ) and [Teixeira &: EmersonI (Il999[ ) added 



Here, the Tycho field stars are plotted as black dots with gray error-bars, the polarization 
targets used by us as black diamonds (those found to be "anomalous sightlines" have been 
over-plotted with X:es), sightlines with r(H20)<0.05 as red open diamonds and those with 
r(H20)>0.1 as filled blue diamonds. No qualitative change in Igo/Iioo vs. A^ is seen at 
the onset of water ice mantles. However, at the optical depth where the 60 and 100 fim 
intensities flatten out as a function of Ay (insert) the linear leo/Iioo vs. Ay relation seems 
to break down, as would be expected if the leo/Iioo ratio was primarily driven by irradiation. 

Together, these results lead us to conclude that we can indeed use the leo/Iioo ratio 
as a tracer of the strength of the radiation field seen by the dust grains. We note that we 
are not attempting to derive physical dust temperatures here, which, due to the complexity 
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of the dust population, would not be reliable from a single ratio. We are simply using the 
FIR ratio as a tracer of the relative strength of the short wavelength radiation in different 
parts of each cloud. In order to minimize any residual effects of variations in the grain size 
distribution, we additionally only use the FIR data to screen the visual extinction data for 
anomalous sightlines, as described in the main text. 
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.04 


119.5±1.2 


HD 110432 


1, 


,40±0, 


.03 


81.0±0, 


.8 


1, 


.59±0, 


.03 


92.1±0, 


.9 


1, 


.70±0, 


.03 


82.6±0, 


.9 


1, 


.73±0, 


.03 


77.3±1, 


.0 


1, 


.66±0, 


.03 


77.2±0.9 


HD 112045 


1, 


.46±0, 


.11 


66.9±3, 


.1 


1, 


.86±0, 


.05 


69.1±1, 


.5 


1, 


.82±0, 


.04 


68.0±1, 


.2 


1 


.87±0, 


.04 


66.8±1, 


.2 


1 


.73±0, 


.04 


68.3±1.1 


HD 112661 


1, 


.49±0, 


.08 


72.4±2, 


.2 


1, 


.80±0, 


.04 


73.7±1, 


.1 


1, 


.86±0, 


.04 


73. Oil, 


.0 


1 


.90±0, 


.03 


70.3±0, 


.9 


1 


.71±0, 


.03 


71.4±0.9 


HD 112954 


2, 


.07±0, 


.06 


41.7±1, 


.7 


2, 


.37±0, 


.04 


42.2±1, 


.0 


2, 


.42±0, 


.03 


40.8±0, 


.9 


2 


.51±0, 


.03 


40.2±0, 


.8 


2 


.29±0, 


.03 


40.2±0.9 


HD 112999 


1, 


.66±0, 


.04 


68.8±1, 


.2 


1. 


.90±0, 


.03 


69.8±0, 


.7 


1. 


.99±0, 


.02 


71.1±0, 


.7 


1 


.95±0, 


.02 


70.1±0, 


.6 


1 


.80±0, 


.03 


70.6±0.8 


HD 113034 


4, 


.03±0, 


.11 


80.0±3, 


.0 


4, 


.64±0, 


.05 


80.8±1, 


.3 


5, 


.05±0, 


.04 


81.7±1, 


.0 


4, 


.97±0, 


.03 


81.8±0, 


.8 


4, 


.49±0, 


.03 


82.7±0.7 


HD 114012 


1, 


.05±0, 


.16 


53.8±4, 


.7 


1, 


.32±0, 


.07 


56.4±2, 


.0 


1, 


.41±0, 


.05 


53. Oil, 


.5 


1 


.53±0, 


.05 


49.1±1, 


.3 


1 


.20±0, 


.05 


54.0±1.5 


HD 114720 


0, 


.83±0, 


.08 


77.7±2, 


.1 


0, 


.88±0, 


.04 


83.3±1. 


.2 


0, 


.97±0, 


.04 


85.5±1, 


.2 





.92±0, 


.04 


79.2±1. 


.2 





.80±0, 


.05 


86.3±1.3 


HD 117111 


1. 


.28±0, 


.03 


76.2±0, 


.9 


1. 


.18±0, 


.03 


77.8±0. 


.9 


1. 


.40±0. 


.03 


77.1±0, 


.8 


1 


.40±0, 


.03 


71.5±0, 


.8 


1 


.38±0, 


.03 


75.5±0.8 



-33- 



Table 3. Coalsack Polarization Curve Fits. 



Star 




Pmaa; 
[%] 




^max 




K 








HD 108417 





68±0 


02 


0.64±0.03 








1.3 




HD 108639 


1 


91±0 


02 


0.52±0.01 





87±0 


13 


0.7 


0.1 


CPD -64 1950 





99±0 


02 


0.63±0.03 








0.8 




HD 108939 


1 


03±0 


02 


0.62±0.05 





42±0 


27 


2.6 


0.2 


HD 109065 





90±0 


01 


0.54±0.01 








10.4 




HD 110245 


1 


22±0 


02 


0.50±0.02 








3.2 




HD 110432 


1 


73±0 


02 


0.62±0.02 





69±0 


12 


4.2 


0.1 


HD 112045 


1 


92±0 


02 


0.57±0.01 








3.1 




HD 112661 


1 


92±0 


02 


0.57±0.01 








0.9 




HD 112954 


2 


52±0 


02 


0.57±0.01 








5.9 




HD 112999 


1 


99±0 


02 


0.56±0.01 





88±0 


10 


2.1 


0.5 


HD 113034 


5 


05±0 


02 


0.574±0.005 








0.35 




HD 114012 


1 


44±0 


03 


0.57±0.02 








2.6 




HD 114720 





95±0 


02 


0.54±0.02 








0.3 




HD 117111 


1 


45±0 


01 


0.58±0.01 








21.6 





Table 4. Stars with data from archives and the hterature. 



Star 


Sp. Class"* 


V 


B-V 


(B-V)o 


V-K 


(V-K)o 




Ay 


ref.*^ Screen? 






[mag.] 


[mag.] 


[mag.] 


[mag.] 


[mag.] 




[mag.] 





Chamaeleon 



Fl 


K4 III (2) 


10, 


,31±0, 


,03 


1 


.56±0, 


.02 


1, 


.39±0, 


F2 


B8 V (2) 


9, 


.90±0, 


.03 


0, 


.51±0, 


.03 


-0, 


.11±0 


F3 


B4 V (1) 


8, 


.03±0, 


.01 


0, 


.42±0, 


.01 


-0, 


.19±0 


F6 


A2 V (1) 


10, 


,85±0, 


,08 


0, 


.74±0, 


.11 


0, 


,05±0 


F7 


B5 V (2) 


10, 


.25±0, 


.04 


0, 


.41±0, 


.06 


-0, 


.17±0 


F9 


KO III (2) 


9, 


.56±0, 


.03 


1, 


.46±0, 


.07 


1, 


.00±0 


Fll 


B9 V (2) 


10, 


.53±0, 


.06 


0, 


.72±0, 


.08 


-0, 


.06±0 


F16 


G2 IV (2) 


11, 


,46±0, 


,03 


1 


.45±0, 


.02 


0, 


,63±0 


F21 


K3 III (2) 


11, 


.39±0, 


.03 


1, 


.94±0, 


.02 


1, 


.27±0, 


F23 


M5 III (3) 


12, 


.93±0, 


.06 


2, 


.55±0, 


.06 


1, 


.63±0 


F24 


B6 IV /V (I) 


7, 


,69±0, 


,01 





.14±0, 


.01 


-0, 


.15±0 


F25 


G8 (3) 


13, 


,23±0, 


,06 


2, 


.22±0, 


.06 


0, 


,94±0 


F27 


M1.5 (2) 


14, 


.80±0, 


.10 


1, 


.26±0, 


.10 


1, 


.47±0, 


F28 


K4(3) 


15, 


.14±0, 


.10 


2, 


.50±0, 


.12 


1, 


.39±0 


F29 


K4 III (2) 


13, 


,33±0, 


,03 


1, 


.74±0, 


.02 


1, 


,39±0 


F30 


K3 III (1) 


11, 


.47±0, 


.03 


1, 


.72±0, 


.02 


1, 


.27±0, 


F32 


FO V (1) 


10, 


.53±0, 


.06 


0, 


.66±0, 


.07 


0, 


.30±0 


F36 


KO (3) 


13, 


.76±0, 


.06 


2, 


.22±0, 


.06 


1, 


.00±0 


F37 


G9 III (1) 


9, 


,58±0, 


,02 


1, 


.27±0, 


.05 


0, 


,97±0 


F39 


K3 III (1) 


10, 


.12±0, 


.04 


1, 


.50±0, 


.15 


1, 


.27±0, 


F40 


B8 III (1) 


7, 


.61±0, 


.01 


0, 


.52±0, 


.01 


-0, 


.11±0 


F41 


B8 V (2) 


8, 


,53±0, 


,01 


0, 


.16±0, 


.02 


-0, 


.11±0 


F42 


A3/4 IV (2) 


8, 


,37±0, 


,01 


0, 


.42±0, 


.02 


0, 


,05±0 


F48 


B9.5 V (1) 


8, 


.88±0, 


.02 


0, 


.12±0, 


.02 


0, 


.00±0, 


F50 


A8 IV (2) 


9, 


,38±0, 


.20 


0, 


.41±0, 


.03 


0, 


.24±0 


F52 


B9.5 V (1) 


8, 


.51±0, 


.01 


0, 


.23±0, 


.02 


0, 


.00±0 



,22 




99+0 


04 


•i 


96+0 


60 

.\J\J 


4 


7+7 


'X 


n 


8+0 


7 


VT! 84 




04 


\ 


.42±0, 


04 


-0 


24+0, 


.12 


Q 


0+n 


, tJ 


X 


8+0 


\ 


VR84 




09 


1 

L i 




09 


-0 




(17 


Q 

o. 




9 


9 


A+O 


1 

, L 


iVl oo 






1 

J_ , 


,UOIIIU, 


OR 


n 


1 zi+n 

. J_4:ZIZU 


07 


9 




A 


1 


.UZHU. 


1 

, J_ 


V xvo^ 


s 


04 


\ 


1 5+0 


O"! 


-0 


.42+0, 


.12 




0+0 

,\J-i—\J. 




\ 


7+0 


\ 


V Xl/O'^ 




16 


4 


54+0 


04 


2 


31 +0 


40 




.3+2, 


2 


2 


5+0 


4 


VR84 




09 


2, 


53+0 


07 


-0 


1 3+0 


.14 


Q 

(J. 


7+0 


R 

.\J 


2 


9+0 


2 


VT! 84 


Q 


05 


4, 


,25+0, 


,04 


1 


.46+0 


.06 


3, 


,7+0, 


,3 


3 


.1+0, 


1 


VR84 




,22 


4, 


.97+0, 


.04 


3, 


.00+0, 


.60 


3, 


.2+1, 


.5 


2, 


.2+0. 




VR84 


g 


,03 


9, 


.09+0, 


.08 


5, 


.96+1, 


.60 


3, 


.7+1, 


.9 


3 


.4+1. 


,8 


W87 


g 


09 


n 


75+0 


03 


-0 

u 


36+0 

.OU— 1— u 


07 


4 


,2+0, 


4 


\ 


.2+0, 


\ 


lV±OkJ 






7 


no+n 

, uynzu. 


, uo 


9 


1 fi+n 

. ±uinu 




A 


9+n 


A 


K 

o. 


4+n 


A 


VV O ( 


s 


iVi 

,Uo 


( . 


.Oz±U, 


• iU 


Q 

o 


.yy±u. 


1 o 
.iz 






d 








\TT> Q A 


s 


33 


8, 


.82+0, 


.10 


3 


.26+0, 


.90 


5, 


.5+2, 


.0 


6, 


.1+1. 


,0 


W87 


s 


,22 


5, 


,67+0, 


,04 


3 


.26+0, 


.60 


7, 


,6+5, 


,1 


2, 


.7+0, 


,7 


VR84 


s 


,11 


4, 


.77+0, 


.04 


3, 


.00+0, 


.30 


4, 


.3+1, 


.3 


1, 


.9+0. 


,3 


VR84 




,03 


2, 


.71+0, 


.07 


0, 


.70+0, 


.06 


6, 


.1+1, 


.3 


2, 


.2+0. 


,1 


VR84 


s 


,24 


7, 


.51+0, 


.08 


2, 


.31+0, 


.70 


4, 


.7+1, 


.1 


5, 


.7+0. 


,8 


W87 


s 


,03 


3, 


,38+0, 


,03 


2, 


.24+0, 


.08 


4, 


,3+0, 


,9 


1, 


.3+0, 


,1 


VR84 




,11 


3, 


.82+0, 


.05 


3, 


.00+0, 


.20 


3, 


.8+3, 


.2 


0, 


.9+0. 


,2 


VR84 


s 


,02 


1, 


.71+0, 


.02 


-0, 


.24+0, 


.09 


3, 


.4+0, 


.2 


2, 


.1+0. 


,1 


MSS 




,04 


0, 


,61+0, 


,03 


0, 


.24+0, 


.18 


1, 


,5+0, 


,8 


0, 


,4+0, 


,2 


MSS 


s 


,07 


1, 


,35+0, 


,02 


0, 


.22+0, 


.16 


3, 


,4+0, 


,8 


1, 


.2+0, 


,2 


MSS 




,02 


0, 


.44+0, 


.04 


-0, 


.07+0, 


.07 


4, 


.6+1, 


.2 


0, 


.6+0, 


,1 


MSS 




,06 


1, 


.08+0, 


.20 


0, 


.57+0, 


.14 


3, 


.2+2, 


.0 


0, 


.6+0. 


,3 


MSS 




,02 


1, 


.01+0, 


.03 


-0, 


.07+0, 


.07 


5, 


.2+0, 


.7 


1, 


.2+0. 


,1 


MSS 





Table 4 — Continued 



Star 


Sp. Class^ 


V 


B-V 


(B-V)o 


V-K 


(V-K)o 


Ry 


Ay 


ref.'' Screen? 






[mag.] 


[mag.] 


[mag.] 


[mag.] 


[mag.] 




[mag.] 





F54 


G6 Ill/TV (1) 


8.40±0, 


.01 


1.09±0, 


.02 


0, 


.89±0, 


.03 




















HD 109753 


B9 ni/iv (1) 

J_J (_/ J- J- J- / J. V \ -L 1 


8.58±0, 


.01 


0.16±0, 


.02 


-0. 


.07±0, 


.05 


HD109885 


B2 HI (1) 


9.02± 0. 


,02 


0.07±0, 


.02 


-0, 


,24±0, 


,06 


HD 110022 


B8 III/IV (1) 


7.84±0, 


,01 


O.lOiO, 


.02 


-0, 


,11±0, 


,05 


HD107875 


B8 V (1) 


9.69±0, 


.02 


0.17±0, 


.03 


-0. 


.11±0, 


.05 


HD 109082 


B9 IV (1) 


8.12±0, 


,01 


0.19±0, 


.02 


-0, 


,07±0, 


,07 


HD109234 


B9.5/A0 V (1) 


9.67±0, 


,02 


0.14± 0. 


,03 


-0, 


,07±0, 


,10 


HD 109565 


B8 H (1) 


9.94±0, 


,03 


0.16±0, 


.03 


-0, 


,11±0, 


,05 


HD 106328 


B4 IV (1) 


8.44±0, 


.01 


0.05±0, 


.02 


-0, 


.19±0, 


,02 


HD109399 


BO. 5 HI (1) 


7.63±0, 


,01 


-0.03±0, 


.02 


-0, 


,28±0, 


,04 


HD 107983 


B9/9.5 V (1) 


8.09±0, 


,01 


0.13±0, 


.02 


-0, 


,07±0, 


,07 


HD 106147 


B9 IV (1) 


8.65±0. 


,02 


0.09±0, 


.02 


-0, 


,07±0, 


,07 


CD- 70 925 


K5 V (1) 


9.16±0, 


.02 


l.lOiO, 


.04 


1, 


.15±0, 


.08 


CPD-69 1677 


AO (1) 


10.06±0, 


,03 


0.27±0, 


.04 


-0, 


,02±0, 


,06 


HD 110118 


B9 IV (1) 


8.44±0, 


,01 


O.lOiO, 


.02 


-0, 


,07±0, 


,07 


CD-69 1024 


FO (1) 


9.75±0. 


.02 


0.44±0, 


.03 


0, 


.30±0, 


.03 


HD 107252 


A3 V (1) 


9.02±0. 


.02 


0.23±0, 


.02 


0, 


.08±0, 


.03 


HD 109233 


AO HI (1) 


9.37±0, 


,02 


0.13±0, 


.02 


-0, 


,02±0, 


,03 


HD 110080 


A5 V (1) 


7.42±0, 


,01 


0.24±0, 


.02 


0, 


.15±0, 


.03 


Ophiuchus 


















001 


G2 V (1) 


8.72±0, 


,02 


0.69±0, 


.03 


0, 


,63±0, 


,03 


005 


A2 (1) 


9.32±0, 


,02 


0.63±0, 


.03 


0, 


.05±0, 


,03 


006 


B9 (1) 


8.60±0. 


.02 


0.39±0, 


.03 


-0. 


.06±0. 


.05 


010 


M4 HI (2) 


13.42±0. 


.02 


2.53±0, 


.03 


1, 


.62±0. 


.02 



2, 


.89±0, 


.02 


2, 


.15±0, 


.05 


4.0±0.7 


0.8±0.1 


MSS 


0, 


.73±0, 


.03 


-0. 


.13±0. 


.13 


4.15±1.16 


95±0 15 

\J • iJKJ —1—\J • 


MSS 


0, 


.36±0, 


.03 


-0, 


,66±0, 


,10 


3 68±0 85 


1.12±0.11 


MSS 


0, 


.54±0, 


.03 


-0, 


,24±0, 


,09 


4.12±1.15 


86±0 10 


MSS 


0, 


.68±0, 


.03 


-0, 


.24±0. 


.09 


3 56± 80 


l.OlibO.lO 


MSS 


0, 


.81±0, 


.03 


-0, 


,13±0, 


.13 


3 95±1 22 


1 04±0 15 


MSS 


0, 


.79±0, 


.03 


-0, 


,13±0, 


,20 


4.81±2.60 


1.01±0.22 


MSS 


0, 


.43±0, 


.03 


-0, 


,24±0, 


,09 


2 75±0 72 


0.74±0.11 


MSS 


0, 


.39±0, 


.03 


-0, 


.49±0, 


.07 


4.00±0.56 


0.96±0.08 


MSS 


-0, 


.01±0, 


.03 


-0, 


,79±0, 


,06 


3 43±0 66 


86±0 07 


MSS 


0, 


.60±0, 


.03 


-0, 


,13±0, 


,13 


4.01±1.61 


80±0 15 


MSS 


0, 


.44±0, 


.03 


-0, 


,13±0, 


,13 


4 00±2 07 


63±0 15 


MSS 


3.07± 0. 


,Uo 


o 

z, 


.oOitU. 


on 


d 




V 


0, 


.67±0, 


.04 


0, 


,00±0, 


,07 


2.60±0.70 


0.74±0.09 


Y 


0, 


.41±0, 


.04 


-0, 


,13±0, 


,13 


3.44±1.66 


0.60±0.15 


MSS 


1, 


.38±0, 


.04 


0, 


.70±0, 


.07 


5.30±1.71 


0.74±0.09 


Y 


0, 


.91±0, 


.02 


0, 


.22±0, 


.08 


5.10±1.37 


0.76±0.09 


MSS 


0, 


.54±0, 


.03 


0, 


,00±0, 


,07 


4.01±1.13 


0.60±0.08 


MSS 


0, 


.76±0, 


.03 


0, 


.38±0, 


.11 


4.58±2.17 


0.42±0.13 


MSS 








1, 


,46±0, 


,03 






VR93 


2, 


.26±0, 


.03 


0, 


.14±0, 


,07 


4.0±0.3 


2.3±0.1 


VR93 


1, 


.47±0, 


.03 


-0. 


.13±0. 


.07 


3.9±0.5 


1.8±0.1 


VR93 


11, 


.18±0, 


.04 


5, 


.10±1, 


.70 


7.3±2.1 


6.7±1.9 


VR93 s 



00 



Table 4 — Continued 



Star 


Sp. Class^ 


V 
[mag.] 






B-V 

[mag.] 




(B-V)o 
[mag.] 




V-K 
[mag.] 




(V-K)o 
[mag.] 








Ay 
[mag.] 




ref.'' 


Screen? 


m A 


V7 TTT n 
r ( 111 [Lj 


12 


.03±0, 


n9 
.uz 


1 
i 


'^9-1-0 
.oZztU 


0"? 

.uo 





.49±0. 


.03 


A 


A 1 -1-0 


0/1 
.U4 


1 


.32±0, 


.11 


A 

4. 


1 -1-0 
.IztU. 


q 
.0 


3 


.4±0, 


1 

. 1 


V nyo 


s 


UiO 


A/T/lo TTT ^^9^ 
iVi^e 111 J 


13 


.59±0, 


09 


9 
Z 


1 -1-0 
.olitU. 


0*? 
.Uo 


1, 


.62±0. 


.02 


O 

O, 


1 Q-l-0 

.lyztu. 


0/1 
.U4 


5, 


.10±1, 


.70 


A 


Q-l-9 
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1 


SM80 




HD 


283855 


A2/0 (2) 


11.40±0 


10 


O.lOiO 


10 


0.05±0.06 


2.21±0 


10 


0.07±0 


14 








2.4±0 


2 


W01/N95 




HD 


283877 


F5/8 V (3) 


9.95±0 


06 


0.57±0 


07 


0.44±0.09 


1.79±0 


07 


1.20±0 


20 


4 


9±4 


5 


0.6±0 


2 


SSS80/N95 




HD 


283879 


B5 V (2) 


11.08±0 


02 


l.OliO 


02 


-0.17±0.06 


2.68±0 


03 


-0.42±0 


12 


2 


9±0 


2 


3.4±0 


1 


WOl 





^Estimated uncertainties, in subclasses, are given in parenthesis 



Spectral cl a sses fr om: KD H94: iKenvon et al.l (19941): M68: iMetrevelil (Il96ah: M S S:lHouk fc CowlevI (119751): N 95 



Slutskii et all (|l980h: SM80: Istraizvs fc MeistasI Jl980[l: SCH85: 'St raizvs et all Jl985[): U8 5: lUngerer et al" 



1984); V R85: IVrba fc Rvdgren (19851 : VR,93: IVrba et al.l (fl993. ): W84: 
200 ll) : YlYal3 (|l997l) 



Whittet et al 



Nesterov et al 



( 1987^ : W87: IWhittet et all dlMZD; WOl: 



(19851): VR84: 



(1995); SSS80: 



Vrba fc Rvdgren 



Whittet et al. 



■^Sightlines that are screened out in the Hipparcos field star based analysis, (see section |4!2|) 
'^The derived value is unphysical and is ignored in the following 

''Based on comparisons of Figure 1 in Arnal et al. ( 1993[) with the Aladin tool ( Bonnarel et al. 2000[) we have made the following target 



identifications in Musca from (H^ZeiaD 11993!) (AMZ): AMZ16=CD-71 836, AMZ18=CD-71 832, AMZ28=CD-70 925, AMZ41=CPD-69 1677, 
AMZ45=CD-69 1024 
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Table 5. Polarization Curve Fits for stars from the literature. 
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0, 
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1 
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.05 


0, 


.570±0. 
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.04 


0, 
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,010 
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0, 
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,020 


0, 
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Musca 
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Table 5 — Continued 



Star 




^max 

[/UJ 






^max 
LA* J 




K 


HD 109753 


2 


12±0 


02 




0.60±0.02 





65±0.15 


HD109885 


3 


82±0 


02 





592±0.005 






HD 110022 


2 


59±0 


01 





587±0.005 






HD107875 


2 


65±0 


03 





505±0.009'^ 






HD 109082 


3 


41±0 


02 





580±0.006 






HD109234 


3 


67±0 


04 





577±0.007 


1 


51±0.18 


HD 109565 


3 


17±0 


04 




0.56±0.01 





93±0.21 


HD 106328 


2 


43±0 


02 





602±0.008 






HD109399 


2 


14±0 


02 





560±0.007 






HD 107983 


2 


36±0 


02 




0.56±0.01 






HD 106147 


1 


83±0 


03 




0.58±0.01 





73±0.18 


CD-70 925 


2 


06±0 


02 




0.54±0.01 






CPD-69 1677 


2 


48±0 


02 




0.59±0.01 






HD 110118 


1 


86±0 


02 




0.58±0.01 





55±0.16 


CD-69 1024 


2 


49±0 


03 




0.58±0.02 





54±0.20 


HD 107252 


2 


42±0 


03 




0.54±0.01 






HD 109233 


2 


10±0 


02 




0.58±0.01 






HD 110080 


1 


42±0 


02 




0.56±0.01 






Ophiuchus 
















005 





50±0 


04 




0.50±0.05 






006 





86±0 


02 




0.65±0.03 






010 


3 


83±0 


14 




0.64±0.04 






014 


1 


42±0 


21 




0.93±0.20 






015 


6 


45±0 


15 




0.89±0.04 






021 


2 


59±0 


07 




0.66±0.03 






024 


4 


31±0 


10 




0.74±0.03 






036 


4 


11±0 


07 




0.77±0.04 






041 


4 


87±0 


13 




0.66±0.04 






042 


3 


65±0 


06 




0.70±0.03 
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Table 5 — Continued 



Star 




[/UJ 










043 


2 


73±0 


07 





71±0 


03 


049 


2 


49±0 


07 





73±0 


02 


050 


2 


34±0 


04 





69±0 


02 


052 


4 


56±0 


14 





71±0 


05 


053 


2 


63±0 


09 





67±0 


05 


055 


5 


53±0 


13 





72±0 


05 


056 


1 


37±0 


41 


1 


05±0 


22 


066 


2 


95±0 


09 





74±0 


06 


069 


1 


96±0 


13 





57±0 


07 


074 


7 


94±0 


17 





72±0 


04 


075 





97±0 


03 





51±0 


04 


080 


3 


99±0 


05 





60±0 


02 


081 


4 


55±0 


10 





50±0 


02 


082 


3 


05±0 


02 





67±0 


01 


084 


5 


47±0 


09 





72±0 


02 


085 


2 


72±0 


06 





76±0 


04 


086 


6 


54±0 


19 





79±0 


04 


087 


4 


26±0 


06 





55±0 


01 


094 


5 


33±0 


09 





69±0 


02 


096 


3 


42±0 


08 





82±0 


06 
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3 


07±0 


04 





66±0 


02 
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5 


97±0 


09 





78±0 


04 


102 


8 


24±0 


34 





66±0 


06 


104 


5 


17±0 


14 





65±0 


03 


105 


4 


86±0 


21 





65±0 


05 


106 


8 


36±0 


12 





61±0 


01 


110 


9 


31±0 


38 





75±0 


06 


111 


9 


07±0 


20 





62±0 


02 


115 


3 


60±0 


25 





55±0 


07 


134 


1 


55±0 


05 





72±0 


05 



K 
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Table 5 — Continued 



Star 


Pmax 


^max 


K 




[%] 


r 1 

[/xmj 




-LOU 


f) Q4+n 09 


f) 71 _|_n no 




144 


X .Ut/ZIZU.U / 


n 61 +n 03 






'4: . O O U . U O 


n c:4+n ni 




147 




65+0 04 

\J . XJ'lJ —i—\J . 






S 40+0 07 


6Q+0 03 

w . vjty _i_w .wo 






7 fiq+n 05 


n fi7+n ni 




R CrA 










1 7q+n 04 


fi7+n 01 




in 


f\7+0 09 


n Qc:+n 03 




12 


f) QO+O 09 


n 77+n 02 




1"^ 


2 77+n 04 


sQ+n 01 

yj .\j>kJ—i—\j .yj ± 




lO 




n 77-i-n ni 




99 


± . ZjtJ _I_U .UU 


n c;n+n no 

U. OU_l_U .UZ/ 






o 1 n_)_n no 

. X W _I_W .\JZj 


n 7s+n ni 






1 8R+n no 


n sn+n ni 




/I o 


1 o 1 _Ln no 

1.81±0.0z 


n '7/1 _i_n n i 

0.74+0.01 




46 


2.61±0.03 


0.81+0.01 




50 


1.12±0.02 


0.76+0.02 




52 


2.01±0.04 


0.70+0.01 




56 


2.09±0.02 


0.67+0.01 




58 


0.84d=0.03 


0.68+0.02 




71 


1.04±0.04 


0.70+0.03 




73 


1.83±0.03 


0.68+0.01 




88 


4.73±0.05 


0.84+0.01 




TY CrA 


0.86±0.02 


0.62+0.04 





Taurus 
HD 28170 
HD 28225 



1.92+0.02 
1.85+0.02 



0.55+0.01 0.96+0.04 
0.56+0.01 0.98+0.05 



Table 5 — Continued 



Star 


Pmai 


^max 


K 




[%] 


[fim] 





HD 


2t 


^975 


3 


16±0 


02 





54±0 


01 





85±0 


04 


HD 


29333 


5 


26±0 


04 





55±0 


01 





88±0 


04 


HD 


29647 


2 


33±0 


03 





74±0 


01 








HD 


29835 


4 


07±0 


03 





50±0 


01 





87±0 


04 


HD 


30168 


4 


07±0 


03 





53±0 


01 





83±0 


03 


HD 


30675 


3 


90±0 


04 





53±0 


01 





95±0 


06 


HD 


279652 


1 


28±0 


05 





48±0 


06 





77±0 


15 


HD 


279658 


2 


82±0 


03 





54±0 


01 





97±0 


05 


HD 


2t 


^3367 


1 


41±0 


02 





70±0 


01 








HD 


2i 


^3637 


2 


72±0 


03 





58±0 


02 





87±0 


09 


HD 


2i 


^3642 


1 


97±0 


01 





61±0 


01 


1 


01±0 


04 


HD 


2i 


^3643 


1 


37±0 


02 





64±0 


01 








HD 


2i 


^3701 


3 


17±0 


04 





60±0 


01 





84±0 


06 


HD 


2i 


^3725 


4 


71±0 


03 





49±0 


01 





71±0 


03 


HD 


2i 


^3757 


2 


89±0 


01 





60±0 


01 








HD 


2i 


^3800 


3 


95±0 


03 





53±0 


01 





85±0 


04 


HD 


2i 


^3809 


6 


69±0 


04 





59±0 


01 





91±0 


03 


HD 


2i 


^3812 


6 


29±0 


03 





55±0 


01 





93±0 


02 


HD 


2i 


^3815 


2 


82±0 


01 





60±0 


01 








HD 


2i 


^3855 


5 


12±0 


04 





52±0 


01 





87±0 


03 


HD 


2i 


^3877 


1 


06±0 


02 





63±0 


01 





96±0 


06 


HD 


2i 


^3879 


4 


14±0 


02 





63±0 


01 









'^For HD 107875 in Musca lArnal et all (119931 ) re- 
port Amaa;=0.641±0.009. This is only case of a sig- 
nificant discrepancy in Xmax for our polarization 
curve fits compared to the literature. This sight- 
line is excluded from the analysis. 
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Table 6. B stars affecting tlie dust 



Cloud 


Star 


Sp. Class 


V 


di 








[pc] 




Chamaeleon 


HD 97300 


B9 V 


9.0 


188±36 


Musca 


HD 109668 


B2 IV 


2.7 


94±4 


Ophiuchus 


a Sco 


Bl HI 


2.9 


225±41 


Opfiiuclius 


p Opli D 


B3/4 V 


6.8 


136±25 


R CrA 


HD 175362 


B3 V 


5.4 


130±16 


Taurus 


HD 30122 


B5 HI 


6.3 


216±34 



Based on Hipparcos trigonometric parallaxes 
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Fig. 1. — Measured multi-band polarization for the Coalsack targets with best-fit polarization 
curve overlaid. Polarization fit parameters are listed in Table |3] 
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Fig. 2. — The wavelength of maximum polarization as a function of visual extinction for the 
newly observed stars behind the Southern Coalsack. The data points have been color coded 
according to the measured values of the total-to-selective extinction along each line of sight. 
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-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 

H-K [mag.] 

Fig. 3. — The J-H vs. H-K color-color diagram for the Coalsack stars shows a consistent 
reddening for all stars (except HD 113034) with a reddening slope around 1.5 (see text). 
The full drawn line shows the intrinsic colors of Main Sequence stars while the dashed line 
is for Giants. 
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Fig. 4. — The measured H-K color excess is plotted (x:es) together with the best fit color 
excesses assuming two reddening laws from the literature. The largest offsets from the color 
excess for reddened intrinsic stellar colors are just above 2a and indicates at most minor 
contamination from circumstellar emission or binary companions. 
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Fig. 5. — The values of Ay and Yiy derived here are compared with thos e extracted 
from the h terature. Data f or Chamaeleon are from IVrba fc RvdgrenI (119841 ). for Ophi- 



uchus from Vrba et al 



Jl993h. for R CrA from I Vrba &: Rydgrenl (119841 ) and for Taurus from 



Whittet et al.l (120011 ). lArnal et aLl (Il993l ) did not pubhsh visual extinctions and their quoted 
Ry values were calculated from Xmax fits, hence Musca is not included here. Good agreement 
is seen for both parameters. 
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Fig. 6. — The wavelength of maximum polarization as a function of visual extinction for 
the sightlines through the Chamaeleon, Musca, Ophiuchus, R CrA and Taurus clouds. For 
Chamaeleon Musca and Taurus a similar structure is seen in the distribution of points as 
that noted for the Coalsack. 
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Fig. 7. — Possible origins of differences between line-of-sight extinction and "effective" ex- 
tinction. For a asymmetrical cloud with its long axis close to the line of sight, the measured 
visual extinction, Ay can both over and under estimate the extinction experienced by the 
gas and dust on the line of sight. See text for details. 
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Fig. 8. — Anomalous sightlines are found comparing the leo/Iioo ratios to the measured visual 
extinction. The top panel for each cloud shows the relationship f o und b etween leo/Iioo ratios 
and Ay for the field star samples selected from the lWright et al.l (120031 ) catalog. The middle 
panels over-plot this best fit on the data for the polarization target samples. The bottom 
panels show the distribution of fit residuals for the field star samples (dashed histograms) and 
polarization target samples (solid histograms). Also plotted are the best-fit Gaussians for 
the field star samples (solid line) and the ±2a limits used for screening targets as anomalous 
(dotted lines). The shaded regions in the top plot for Chamaeleon illustrate the approximate 
areas corresponding to the "sightline types" discussed in the text and illustrated in Figure 
[3 This "debiasing" technique fails for the Coalsack sightlines, presumably due to the strong 
FIR background in the Galactic plane. See text for details. 
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Fig. 9. — For the Coalsack, establishing a nominal Ay vs. leo/Iioo slope fails, presumably 
due to the influence of background emission in the IRAS data. We therefore screened the 
Coalsack data based on similarities with the Chamaeleon and Taurus plots. The screened-out 
sources are marked with X-es. 
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Fig. 10. — The offset from the best-fit leo/Iioo vs. Ay relations (Tycho selected field stars) 
are plotted as functions of the 2- dimensional distance from nearby hot stars. For lines of 
sight through Chamaeleon the dominant star is HD 97300 (top left); for Musca HD 109668; 
for Ophiuchus, a Sco and p Oph D both yield noticeable effects. For lines of sight through 
R CrA and Taurus the dominant stars are HD 175362 and HD 30122, respectively. Stars 
from the Tycho field star samples are plotted with x:es while diamonds indicate stars in 
our polarization samples. Open diamonds indicate sightlines classified as "anomalous". The 
polarization sample stars with the largest offsets tend to be located close to the hot stars, 
indicating that many of the "anomalous" leo/Iioo ratios are due to the additional irradiation 
from these stars. 
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Fig. 11. — Plots of Xmax VS. Ay for the debiased sample for Chamaeleon, Coalsack and 
Musca. For each cloud, the best fit (using a robust linear fit) is over plotted as a full drawn 
line. The dashed lines represent the "universal slope" of 0.028 and zero intercepts determined 
from the average Ry for the cloud. The Hipparcos screened sample yield similar plots 
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8 10 



Fig. 12. — Same as Figure[II]for Ophiuchus, R CrA and Taurus. For Ophiuchus and Taurus 
the Hipparcos screened samples yield similar plots. For R CrA a significantly different 
selection is encountered which is illustrated in Figure | 
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Fig. 13. — Plot of \max VS. Ay for the debiased R CrA sample screened based on the 
Hipparcos field stars. As in Figure [TT], the best fit (using a robust linear fit) is over plotted 
as a full drawn line. The dashed line represents the "universal slope" of 0.038 and zero 
intercept determined from the average Ry for the cloud. 
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Fig. 14. — The slope (left) and zero intercepts of the Xmax vs. Ay (right) for the six 
clouds under study are compared to the average of the total-to-selective extinction, using 
the Tycho screened samples. The full drawn and dashed lines in the left-hand panel illustrate 
the weighted average of the slope and its uncertainty. The full drawn line (right-hand panel) 
shows the best fit for Xmax vs. Ay, allowing a non-zero intercept. The short-dashed line 
shows the best fit assuming a zero intercept. The long-da shed line corresponds to the Ry 
vs. Xmax correlation found by lWhittet fc van Bredal (119781 ). 
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Fig. 15. — The adjusted location of the peak of the polarization curve is shown as functions 
of alignment efficiency (pmax/Ay; left panel) and Ry (right panel) for the Tycho screened 
sample. In both cases the measured \max has been adjusted to the Chamaeleon value by sub- 
tracting the difference in derived \max{,-^v = 0) between the cloud and that for Chamaeleon 
for the sightlines in each cloud. I.e. A;,„,=A^,,-(A„„,(Ay=0)'='°"^-A„„,(Ay=0)^'^'^'"). 
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Fig. 16. — The offsets from tlie Xmax vs. Ay fits are plotted as a function of tlie 1/0 H2 
temperatures for six stars in Opfiiuclius (filled diamonds) and one each in Chamaeleon (open 
gray diamond) and the Coalsack (open black diamond). No correlation is evident. 
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Fig. 17. — The offsets from the fits in Xmax vs. Ay are shown in histogram form (full drawn 
histograms and hnes for the Tycho screening, dashed histograms and lines for the Hipparcos 
screening). The offsets are from the "universal slope" relations. 
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Fig. 18. — The scatter in the fit of Xmax vs. Ay is plotted against the bright source fraction 
of the embedded young stellar objects (Open symbols for the Hipparcos screened sample and 
filled diamonds for the Tycho screened sample). The curves are the best fit exponentials. 
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Fig. 19. — The combined binned fractional bolometric l uminosity function of YSOs in 
Chamaeleon Ophiuchus and Taurus is shown as diamonds (jChen et al.l 119971 ). This shows 
the fraction of the total YSO luminosity, per unit luminosity, originating from a given sub 
sample of YSOs. The binned fractional X-ray luminosity function of Ophiuchus YSOs is 
plotted as histograms. The binned fractional luminosity functions peak at ~l-2 
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Fig. 20. — The ratios of I12/I100 (blue), I25/I100 (green) and Igo/Iioo (red) are plotted for 
Chamaeleon, Musca, Ophiuchus, R CrA and Taurus, against the on-the-sky distance to the 
dominant nearby hot star. As discussed in the main part of the text, the 60 fim to 100/im 
ratio responds noticeably to the presence of the hot, high luminosity star. The 25 to 100 
fim ratio shows a smaller response, while the 12 to 100 fim ratio shows very little response 
to the presence of the hot stars. This is consistent with an origin in dust heating through 
irradiation of the grains, but would be difficult to explain in terms of changes in the relative 
abundances of the different grain populations. 
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Fig. 21. — The ratio Im/ Imn vs. Ay is plot t ed fo r Taurus, with the sighthnes from 
Murakawa et al.l (120001 ) and iTeixeira &: EmersonI (119991 ) added. Black dots with grey error- 
bars represent the Tycho field stars, black diamonds the polarization targets use in this paper 
(sightlines screened out as being anomalous are indicated with an X), red open diamonds 
represent sightlines with measured r(H2O)<0.05 and filed blue diamonds r(H2O)>0.1. The 
insert shows the 60 and lOO/im intensities for each sightline. 
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Fig. 22. — The dispersions around the best leo/Iioo fits (Figure [H]) using the Tycho based 
screening are plotted as a function of the bright star fraction for each cloud (we have also 
added the results for a similar analysis of the Lupus I cloud). A very well defined sequence 
is found, well fit by a power law function. 



